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Although the two-body charmed decays B( s ) — > and D^V, where P(V) denotes a light 

pseudoscalar (vector) meson, are CKM suppressed comparing with the B( s ) — > _D|*jP and D^V 
decays, they are important in the CKM angle 7 = (^3 extraction method. We investigated these 
decays in the perturbative QCD approach to the leading order of mu/ms and Aqcd/ttid expansion. 
We find that the nonfactorizable emission diagrams and the annihilation diagrams are not negligible 
in many of these channels. The numerical results show that most channels have branching ratios 
with an order of 10 -6 or 10 -7 . The ratio needed for the CKM angle 7 extraction is estimated 
as r = M(f^5njpjj = 0- 092 io'.oo3-o'.oo3 - which is too small for the experiments. Some of the 
-B( s ) — > D( S )V decays have a very large transversely polarized contribution that can reach 80%. 



PACS: 13.25.Hw 



I. INTRODUCTION 

The study of B physics plays an important role in precise test of the standard model, extraction of the Cabbibo- 
Kobayashi-Maskawa(CKM) matrix elements, searching for the origin of CP violation and new physics signals and even 
discovery of new hadronic states. The hadronic B decays offer an opportunity to understand the nonperturbative 
QCD. After years of hard work, theorists develop many approaches to deal with the nonleptonic decays of B mesons, 
such as QCD factorization approach soft collinear effective theory Q], perturbative QCD approach(PQCD)Q], 
QCD sum rules[4|], light cone sum rules[5j. At the experimental side, the two B factories have accumulated a great 
amount of data, which can be used to either test various theoretical approaches or carry on new physics investigations. 
The two body charmed decays of B mesons B^ — > D K ,JP and .D^V^ which are important in the extraction of the 



CKM angles [6[, have been investigated in the PQCD approach [7|, |8|, |9(. These channels are induced by the b 



transitions, which are CKM favored. However, the method of CKM angle 7 extraction [lfj also requires another 
category of charmed meson B decays, which are induced by b — > u transition. The interference between the b — > c and 
b — > u transitions gives the measurement of the CKM angle 7. These b — * u decays are CKM |V^i,/Vc&| suppressed, 
thus will have smaller branching ratios. In this paper, we will investigate these B(B S ) — > D^P, D^V, D^P, D*^V 
decays in the PQCD approach. 

Unlike the collinear factorization in the QCD factorization approach and soft-collinear effective theory, the kx 
factorization is utilized in the PQCD approach. In this approach, the transverse momentum of valence quarks in the 
mesons are kept to avoid the endpoint singularity. Therefore, only in this factorization method, one can calculate the 
color suppressed channels as well as the color allowed channels in charmed B decays. The conventional non-calculable 
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annihilation type decays are also calculable in the PQCD approach, which is proved to be the dominant strong phase 
in B decays for the direct CP asymmetry In the PQCD approach, the most important uncertainties come from 



the hadronic wave functions. We will use the same hadronic parameters determined from the charmed B decays 
induced by b — > c transitions {J to reduce uncertainties. The numerical results show that the b — > u transition decays 
are indeed heavily suppressed comparing with those b — > c transitions. Thus the CKM angle measurement method 
do face difficulty experimentally 

The paper is organized as follows: In Scc.|TT]we list the analytic formulae needed in the calculation, including the 
Hamiltonian and definition of the momenta, PQCD factorization formulae of all the diagrams, and the expressions 
of the amplitudes for all the decay channels. Sec. IIIII contains the numerical results we obtain, and some discussions. 
A brief summary is given in Sec. IIVI The wave functions, decay constants and some functions that appear in PQCD 
approach are put in the appendices. 



II. ANALYTIC FORMULAE FOR DECAY AMPLITUDES 

In B meson weak decays, there are three natural energy scales involved: W boson mass, b quark mass scale and 
the hadronic scale Aqcd- The electroweak physics higher than W boson mass can be perturbatively calculated. The 
physics between b quark mass and W boson mass can also be calculated using the rcnormalization group equation. 
Both of these two contributions are included in the well calculated effective Hamiltonian of the four quark operators. 
The physics below b mass scale is complicated, where we have to utilize the factorization theorem to factorize the 
non-perturbative contribution out. In this purpose, we do the l/m^ expansion in the soft collinear effective theory 
Q. Unfortunately, there are some contributions, which produce endpoint singularity. In order to deal with this 
singularity, usually there are two ways of doing it. One of them is the PQCD approach, in which we keep the 
transverse momentum of the valence quark. By doing this kx factorization, a new series of double logs are generated. 
Using the renormalization group equation, we resum these logs to give a Sudakov factor, which suppresses the endpoint 
contribution. Finally the decay amplitude becomes 

M= f d 4 fc 1 d 4 fc 2 d 4 fc30B(fci,t)TH(fc 1 ,fc 2 ,fc 3 ,t)0p 2 (fc 2 ,t)0p 3 (fc 3 ,i)e s(fe -*\ (1) 



where Th is the hard part that is perturbatively calculable, and 4>M a- re the hadronic meson wave functions that 
is non-perturbative. The Sudakov factor e s ( k i^) resulting from the resummation of double logarithm, relate the 
perturbative and non-perturbative dynamics. 

In charmed B decays, we have one more intermediate energy scale, the D meson mass. Therefore, one encounters 
another expansion scries of mu/wB. The factorization is only approved at the leading of to^/tob expansion [zlS], 
so that we will keep only the leading order contribution numerically, unless explicitly mentioned. 

A. Effective Hamiltonian and Kinematics 

For the processes considered, only tree operators contribute, and the effective Hamiltonian is given by 

H eff = ^Vu b V L * D [CMOM + CaMOaOO] , (2) 
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with D — d, s, and 



Oi = {u a bp)v-A{Df}C a ) V -A ,02 = {u a ba)v-A{DpCf3) V -A ■ (3) 

Here, a (3 are the color indices, (q~iq 2 )v-A = <Zi7 M (l — 7 5 )<Z2- And the V u bV* D are the corresponding CKM matrix 
elements. 

Sandwiching the above Hamiltonian between the initial and final state mesons and factorizing the matrix elements, 
the combinations of the Wilson coefficients usually appear. Conventionally, they are defined as 

a t = 2 + 0i/3 , a 2 = 0! + 2 /3, (4) 

where, a\ and a 2 correspond to the color favored and color suppressed contribution, respectively. The light-cone 
coordinates are used in this paper, with which the vector V is expressed as V — ( — , v , Vj_), where Vj_ = 
(V 1 , V 2 ). The momenta of B, D and the light mesons are respectively Pi, P 2 and P3, which are defined as 

P 1 = ^|(1,1,(U), P 2 = ^(l,r 2 ,(U), P 3 = ^|(0,l-r 2 ,O i ), (5) 

with r = rriu/mB- The momenta of the light quarks in B and D mesons are denoted by k\ and k 2 , respectively, 
whereas k% represents the momentum of the quark in the light meson. Both k^ and k^[ contribute, but considering 
k 2 ~ 0(A) and ki ~ 0(o;itob), we drop the term k\ ■ k\ for the emission diagrams and k\ ■ k 2 for the annihilation 
diagrams. The effect is equal to dropping k^ for emission diagrams and dropping k\ for annihilation diagrams. Their 
explicit expressions are 

k\ = {x\ — 0, kijj for emission diagrams, 
k\ = (0,xi — ^-,kuj for annihilation diagrams , 

k 2 = ( a!a ^,0 > k 2± ) , fc 3 - (0,x 3 (1 "^ )mB ,k 3± ) . (6) 
Here, x\, x 2 and 23 are the momentum fractions, and ku_, k 2 j_ and k 3 j^ are the transverse momenta of the quarks. 



B. Factorization formulae of B — > DP 

The contributions to the B — > DP processes can be divided into two types: The emission diagrams(scc Fig. [1] 
denoted by a subscript int in the following formulae), in which the light quark in the B mesons is hadronized into 
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FIG. 2: Annihilation diagrams in pQCD approach for B — * DP decays. 



one of the final state mesons as a spectator, and the annihilation diagrams (sec Fig. [5J denoted by a subscript exc 
in the following formulae) with no spectator quarks. The first two diagrams of Fig. [pare the factorizable diagrams, 
which can be factorized as a product of the decay constant of D meson and a B to light meson transition form factor. 
Summing the expressions of these two diagrams together, we obtain 

rl/A 



6nt(ai) = 8nC F fD I dxidxz 

lo Jo 



bidbib 3 db 3 4> B (xi,bi) 



x{[(2 - x 3 )Mxs) - r (l - 2x 3 ){^ P {x 3 ) - <fe(x 3 ))] 
xaiit^Eiit^hiixt, (1 - x 3 )(l - r 2 )MM)St{x 3 ) 



+2r <f p (x 3 )a l {tf ) )E l {t [2} )h l {l - x 3 , Xl (l - r 2 ), b 3 , h)S t ( Xl )} , 



,(2) 



(7) 



where ai is the corresponding combination of Wilson coefficients and ro = mo / tub , with mo as the chiral mass of 
the pseudoscalar mesons. The expressions of the PQCD factorization functions hj, jet function St(x) and Ej(t l j) and 
scales t l j, with j = i, a, d, / and I = 1 or 2, are listed in Appendix [Bl 

In the factorized diagrams of annihilation contributions (the first two diagrams of Fig. [J), the B meson is factorized 
out. And the combination of these two diagrams give 

,1 ,1/A 

Cexc(ai) = 8irC F f B / dx 2 dx 3 / b 2 db 2 b 3 db 3 (p D (x 2 ,b 2 ) 



x (-x 2 0p(x 3 )-2ror(l+s 2 )^(x 3 ))a l (tW)i; a (ti 1 ))/ la (l-a;3,a; 2 (l--r 2 ),63,6 2 )5 t (x3) 
+ ((1 - x 3 )0 P (x 3 ) + r r((2x 3 - 1)0? (x 3 ) + (3 - 2x 3 ) ( / )p (x 3 )))a l (t^)E a (t^) 
xh a (x 2 ,(l-x 3 )(l-r 2 ),b 2 ,b 3 )S t (x 2 )] . (8) 

The last two diagrams of Fig. [1] and [3] arc nonfactorizable diagrams. Generally, if the two final state mesons are 
both light ones, the nonfactorizable contributions of the emission diagrams are very small, because of the cancelation 
between the two nonfactorizable diagrams. While for D mesons, since the heavy c quark and the light quark is not 
symmetric, the nonfactorizable emission diagrams give remarkable contributions. The expression of the nonfactorizable 
emission contributions is 

l/A 



(9) 



A^i n t(ai) = 16tt^2N c C f / [dx] I bidbib 2 db 2 (j)B{xi,bi)(j) D (x 2 ,b 2 ) 

Jo Jo 

x Oe 2 <M.t 3 ) + n>(x 3 - 1)(#,(X3) + 0p(a*))) ai(t<P)E d (t<p)h$\xi,bi) 

+ ((x 3 + x 2 - 2)<f> P (x 3 ) + r (l - x 3 )(r p (x 3 ) - <jPp{x 3 ))) Oi^E&f^f&iM) 



and that of the nonfactorizablc annihilation contributions is 

/-l /-i/A 

M exc {ai) = 16ny/2N c C F / [dx] / bidbib 2 db 2 4> B (xi, b\)(j)r>{x 2 , b 2 ) 

Jo Jo 

x [((x 3 - 1)0p(.t 3 ) - r r((a: 2 - x 3 + 3)<$,(> 3 ) + {x 2 + x 3 - l)<j>l{x 3 )))ai{tf)E f {tf)hf{x 
+ Mp{x 3 )-r Q r{{x 3 -x 2 - l)<j> P {x 3 ) + (x 2 + x 3 - l)tt>%(x 3 )))c^(tf ) )E f (tf ) )hf\x i ,b i )] ■ 
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C. Amplitudes for B — > DP 

With the functions obtained in the above, the amplitudes of 22 B — ► Z)P decay channels can be given by 

A(B-^D\-) = ^m%V ub V c * d {^nt{a2)+M mt (C 2 /S)+^ xc {a 1 )+M exc {C 1 /3)) 1 (11) 

A{B-^D°K-) = ^m%V ub V c * s (C m t(a 2 )+M l nt{C 2 /3) + ^ xc (a 1 )+Me X c(C 1 /3)), (12) 

A(B- -> ^-7r°) = ^Lm%V ub V* d ^= (£ in t(ai) + M^d/S) - UM - M exc (Ci/3)) , (13) 

.4(£T ^ D-/?°) = ^L m %V ub V c * s {£ eX c(ai)+M eX c(Ci/3)), (14) 
v2 

^(B" -> D~7] nfi ) = ^m%V ub V c * d ^= (6nt(ai)+Mnt(Ci/3)+^xc(ai)+M e xc(Ci/3)), (15) 

^(B-^^ s -7r ) = ^m^KtVS tent(ai)+M int (Ci/3)), (16) 

^l(B-^P s -A" ) = ^m%V ub V^(^xc(ai) + Mexc{Ci/3)), (17) 
v2 

.4(5" -» D~rj n n) = -^m%V ub V* s -^= (&„ t (ai) + M mt (Ci/3)) , (18) 

.A(-B _ — * Djr) S s) = ^m%V ub V c * s ^e X c(ai) + M exc (C 1 /3)), (19) 

.4(5° -> 5°7T ) = ^| TO |y u6 v;* d -^ (-&„ t (a 2 ) - M mt (C 2 /3) + Ce, c (a 2 ) + X e:cc (C 2 /3)) , (20) 

„4(5 ^£°A ) = ^m%V ub V c * s {&nt(a 2 ) + M in t(C 2 /3)), (21) 

.4(5° D°n nn ) = ^Z m %v ub V c * d ^= (Ut (02) + MintiCi/d,) + ^ exc (a 2 ) + M exc (C 2 /3)) , (22) 

^(B°^Z?-7r+) = ^m|KfeKdtent(ai)+M i „ t (C' 1 /3)+^ c (a2)+^exc(C2/3)) 1 (23) 

„4(S° ^ .D S 77T+) = ^|m|K6K* s tenf(ai)+Mnt(Ci/3)), (24) 

A(B°^D-K+) = ^m%V ub V c * d ^ exc (a2)+M exc (C 2 /3)), (25) 

.4(5° -> £ 7r ) = ^ m |K fe y c * -L (^ xc (a 2 ) + M e *c(C 2 /3)) , (26) 

„4(5 S ° ^ 5°A'°) - ^m|KbKM&nt(a 2 )+M m i(C 2 /3)), (27) 

^(Bj-^DVfi) = ^m%V ub V c * s ^={Uc{a 2 )+M exc (C 2 /3)), (28) 

,4(5° ^Vs) - -^mlKftKs (6nt(aa) + -M^t (C 2 /3)) , (29) 
v2 

^(S s °^^-7r+) = ^m^KiV^ (6^(02) +AW<V3)), (30) 

A(B° -> D-K+) = ( ^m A B V ub V: d ^ mt {a 1 )+M lnt {C 1 /3)), (31) 

A{B° S -> £>7A'+) = -^m%V ub V* s {Cint(ai) + M inf (Ci/3) + + M exc (C 2 /3)) , (32) 
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where r) n n = -^={uu + dd) and rj sS = ss. We treat r\ and rf as the mixtures of r\ nfl and rj sS with 



where 6 = 39.3° ± 1.0 C 






(33) 



D. Decay Amplitudes for B — > D*P and B — > £r*'V modes 



Due to the conservation of the angular momentum, only the longitudinally polarized vector mesons are generated 
in the B — > I?V and B —> D*P decay modes. According to the similarity between the Lorentz structures of the wave 
functions for the pseudoscalar mesons and vector mesons, the emission diagrams' factorization formulae of B —> DV 
modes can be obtained by making the following substitutions in Eq. {7} and (J9]) up to the leading order and leading 
power of 1/mp: 



^, 4>l 



,1 ,1/A 

SnC F fB / dx 2 dx 3 / b 2 db 2 b 3 db 3 4> D (x 2 , b 2 ) 
Jo Jo 



-<py, r -> ry, /p -> jV- 
The factorization formula for the annihilation diagrams of i? — > decays are listed below: 

Cexc(^z) 

x {-x 2 <j> v {x 3 ) + 2r v r{l + x 2 )cf> v (x 3 )) ai (t^)E a (t^)h a (l - x 3 , x 2 {\ - r 2 ), b 3 , b 2 )S t {x 3 ) 

+ ((1 - x 3 )4> v {x 3 ) + rvr((l - 2x 3 )<^(x 3 ) - (3 - 2x 3 )<^3))K(^ 2) )£a(^ 2) ) 

x/i a (x 2 , (1 - x 3 )(l - r 2 ), 6 2 , 6 3 )S t (x 2 )] . 
/■! /-i/A 

167tV2AcC_f / [efx] / b 1 dbib 2 db 2 (j)B{xi,bi)4>D(x 2 ,b2) 
Jo Jo 

x[((x 3 - l)0v(a:3) + rvr((x 2 + x 3 - l)<fc(x 3 ) + (x 2 - x 3 + 3)<^(x3)))ai(tf ) )Ef(tf ) )hf\x i , h) 



(34) 



(35) 



M exc (ai) 



+ (x 2 0y(x 3 ) - r v r{{x 2 - x 3 + l)<j> v {x 3 ) - {x 2 + x 3 - l)<j> v {x 3 )))adtf)E f {t\ z, )h) z \x l ,b i )] ■ (36) 

By changing the pseudoscalar nonet to corresponding vector nonet in Eq. (|ll|) -(|32 |) . one will obtain the decay ampli- 
tudes formula for the corresponding B — > DV mode. 

While in order to obtain the factorization formulae for B — > D*P mode, the substitutions 

§d -> </>p, fo -> /r>*, "id -> "ir>* (37) 
should be made in Eq.((7|) and j9j for emission diagrams. The annihilation diagrams' formulae are listed below: 

£exc(&i) 

x (-x 2 0p(x 3 ) + 2r o r(l-x 2 )^(x3))a l (t«)^ a (ii 1) )/ la (l-X3,x 2 (l-r 2 ),63,6 2 )5 t (x 3 ) 

+ ((1 - x 3 )0 P (x 3 ) + r r(<^(x 3 ) + <^,(x 3 ))) ai (tf)E (4 2 ))/i (a; 2 , (1 - x 3 )(l - r 2 ),b 2l b 3 )S t (x 2 )] (.38) 

A'(exc(ai) = \&-Ky/2N c CF / [dx] / bidbib 2 db 2 4>B(xi,bi)4>D(x 2 ,b 2 ) 



*( 2 )\f,( 2 ) ( 



fl /-1/A 

87tCf/s / dx 2 dx 3 / b 2 db 2 b 3 db 3 (j)D(x 2 ,b2) 



n 



x [((x 3 - l)0 P (x 3 ) - r r((l - x 2 - x 3 )0 p (x 3 ) + (1 - x 2 + ais)^^)))^^)^/^)^^,^) 



+ (x 2 P (x 3 ) +r r((x 3 -X2- l)<£p(x 3 ) + (x 2 + x 3 - l)^,^)))^^)^/^)^^.^)] 



t ( 2 hf,( 2 ), 



(39) 
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The decay amplitudes for these modes can be got through changing the D mesons to corresponding D* mesons in 

lHU-(32). 

The situation for B — * D*V mode is a little more complicated. Both the longitudinal polarization and the transverse 
polarization contribute. The longitudinally polarized factorization formulae are obtained by making the following 
substitutions in Eq. ([7|) and ((9]): 

<I>p -> -<f> v , 4> P -> 0vj 0V"j r o -> rv, /p -> /v, 0n -> /d -> /d*, -> m D .. (40) 

The annihilation diagrams' formula are 

/•l /-l/A 
^oxc(ai) = SnC F f B / dx 2 dx 3 / b 2 db2b 3 db3<f> D (x2,b 2 ) 
Jo Jo 

x (-x 2 y (a;3) + 2ryr(x 2 -l)^(x 3 ))a J (^ 1) )S a (ti 1) )/i a (l-a;3,x 2 (l-r 2 ),63,&2)^(x3) 

+ ((1 - x 3 )0y(x 3 ) - tvr(0V(s 3 ) + ^(^3)))a i (4 2) )K(ii 2) )/ia(^2, (1 - z 3 )(l - r 2 ), &2, & 3 )St(x 2 )] (41) 

/•! /-l/A 

M cxc (ai) = -16tt\/2N c Cf / [da;] / b 1 dbib2db 2 (p B {xi,bi)^D(x2,b 2 ) 

Jo Jo 

X[((x 3 ~ 1)<M^ 3 ) - ryr((x 2 -x 3 - l)^(x 3 ) + (ar 2 +^3 - l)4fy(x 3 )))c*(tf ) )E f (tf ) )hf\x i ,b i ) 
+ (x2<t>v(x 3 ) + r v r((l -X2- x 3 )0 s v (x 3 ) + (x 2 - x 3 + ^{x^a^tf^^tf^f^h)] ■ (42) 

The longitudinally polarized decay amplitudes A N s can be obtained through replacing the pseudo-scalar mesons by 
the corresponding light vector mesons and the D* meson in (|ll [l - (|32j) . 

The transversely polarized amplitude of B — > D*V mode can be decomposed as 



A T (el*,el*) = ^m B V ub V c * D 



(43) 



where A s/p = £,llt + ^Int + &Lc + M s J x P c- , M^, Cxc, Mlic correspond to the factorizablc emission diagrams, 
nonfactorizablc diagrams, factorizablc annihilation diagrams and nonfactorizablc diagrams respectively, and their 
analytic expressions are given in Eq. (|44p -Eq. (|5ip . and ty are the respective transverse polarization vectors of D* 
and vector mesons, n is the light cone vector in which direction the momentum of the vector meson is defined and n 
is the opposite direction. e^ pa is the antisymmetric tensor with g 0123 = 1. 

The transversely polarized contributions arc suppressed by r or rv, and their expressions are given by 

,1 ,1/A 
£int( fl i) = %nC F f D r \ dx Y dx 3 I b 1 db 1 b 3 db 3 (j) B (x 1 ,b 1 ) 
Jo Jo 

x{[-0y(z 3 ) + r v ((x 3 - l)^(x 3 ) + (x 3 - 3)4> v v )} 
xaiit^Eiit^hixx, (1 - x 3 )(l - r 2 ), b 1 , b 3 )S t (x 3 ) 
+rv[<P v (x 3 ) - ^(X3)]a i (4 2) )^(4 2) )^( 1 - ^,xi(l - r 2 ), b 3 , &i)S t (zi)} , (44) 

,1 ,1/A 

^Sit(°») = %nC F f D r \ dxidx 3 / bidbib 3 db 3 4> B {xi,bi) 
Jo Jo 

x{h0y(^ 3 ) ~ r v {{x 3 - lW v (x 3 ) + {x 3 ~ 3)^)] 
xoi^)^^ 1 *)^, (1 - x 3 )(l - r 2 ), h, b 3 )S t (x 3 ) 
+r V W v (x 3 ) - cj ) l(x 3 )]a l {tf ) )E l {tf ) )h{l - x 3 , Xl (l ~ r 2 ),b 3 , b 1 )S t (x 1 )} , (45) 



1/A 

16TTy/2N c C F r I [dx] / hdbib 2 db 2 (f>B(xi,bi)4%(x2,b2) 



-x 2( t>l{xMt^)E d {t^)h { p{ Xl A) 



+ (-rv(<ft(xa) + (^2 + 2x 3 - 5)^(x 3 )) + (x 2 - 3)<fc(x 3 )) ai {tf)E d {tf)hf (x,A) , (46) 



fl /-l/A 

167r\/2iV c C_Fr / [da;] / bidbib 2 db 2 (j)B{xi,bi)(f)^ ) {x 2 ,b 2 ) 



o 



£cxcK) 



-x 2 ^(x 3 )a l (4 1) )^(4 1) )^ 1) ( 2 : 4 ,6 l ) 
+ (rv(<^(x 3 ) + (2x 2 + 2x 3 - 5)^(x 3 )) + (x 2 - 3)^(x 3 )) ai (tf)E d {tf)hf \ Xi M) 



rl /-i/A 
%irC F f B I dx 2 dx 3 / & 2 d& 2 & 3 d& 3 0^(x 2 , 6 2 ) 
Jo Jo 

rry f(x 2 - l)^(x 3 ) + (*2 + l)^(x 3 )Vi(*a ^^(^^^(l - z 3 ,x 2 (l - r 2 ), b 3 , b 2 )S t (x 3 ) 



(47) 



+ (jr y (x 3 (<^(x 3 ) - ^(x 3 )) ~ 2<^(x 3 )) + r 2 ^(x 3 ) 
xa,:(4 2 ))K(tl 2) )/i n (^, (1 - z 3 )(l - r 2 ), 6 2 , & 3 )S { (x 2 ) 



(48) 



/•l /-l/A 

87rC F /ij / dx 2 dx 3 / 
Jo Jo 



b 2 db 2 b 3 db 3 (p D (x 2 ,b 2 ) 

-rr v ^ 2 - lW v {x 3 ) + (x 2 + l^fa^Oiit^Eait^haQ. - x 3 , x 2 (l - r 2 ), b 3 , b 2 )S t {x 3 ) 
+ (rr v (x 3 (^ v (x 3 ) - <^(x 3 )) + 2<j> a v (x 3 )) + r 2 <$(x 3 ) 
x ai {t^)E a {t^)h a (x 2l (l- X3 ){l~r 2 )MM)St{x 2 )\ , (49) 



bidbib 2 db 2 (/)B(xi,b 1 )4)] D (x 2 ,b 2 ) 



A^ xc (a 4 ) = lfo^/2N c C F / [dx] / 

Jo Jo 

x^rrv^lix^aiit^Ef^hfixiM) 
H-X2r 2 <Pl{x 3 ))a t {tf)E f {tf)hf (x^k)] . 

i fl /-l/A 

ML(a<) = 16^V2N C C F / [dx] / M&i&adMflOci.&OtfSte.&a) 

Jo Jo 

x[(-2rr y ^(x 3 )a l (4 1 ))£; / (4 1) )4 1) (^,6 4 ) 
+(-x 2 r 2 ^(x 3 ))a i (4 2) ) J B / (4 2) )4 2) (x i ,6 i )] . 



(50) 



(51) 



III. NUMERICAL RESULTS AND DISCUSSIONS 



Although the meson wave functions are not perturbatively calculable, they are universal for all the decay channels. 
We can determine them from the well measured decay channels, such as the B — ► D^P decays Q]. We use the 
same light cone distribution amplitudes as those wc obtained in Ref. Q]. The Lorentz structure of wave functions, 
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TABLE I: Branching ratios of B/ s \ — ► DP decays calculated in PQCD approach together with experimental data[3| (unit: 

HT 6 ). 



Modes class Experiments Our results 



B° 


-> D77T+ 


T 


15.3 ± 3.5 


t-» „-L1 7 A -LI Q4-1 9 
QQ ri-|--L / .y-|-o. z 

OO.U_ 12 .6_i.g_3.i 




B° 




E 




-0.15-0.12-0.04 J x 1U 


_ 2 


B° 


fSO 
— > D 7T 


C 




(4.50±?J±i;l±°.lt) x 10 


■> 


— n 

B" 


-» D K 


c 




- — ^-un Q~_i_n 7R-i-n is 
L.ty_ Q 74-0.55-0. 17 




B° 




c 




z _ _,,_|_~1_|_1 "_Lfl R1 \ . _ 

(6.i6±i:|±i:|±g-.s^) x io 


2 


B u 


-> D u rj 


c 




/ . i d+2.1+1.0+0.41 \ 1 a- 
( 4 -^- ; -l'5-0.90-0 .38 ) X - U 


. 2 


B° 


-> D~-k + 


T 




_ „_i_n _ w_i_n (17-i-n 1 1 
±.±U_q 43_o.07-0. 10 




B~ 


— n 


c 




_ _ ___j_n n7j_n n _ j_n ri9 
u --' -0.06-0.03-0.02 




B" 


-> D ^" 


c 




-• oy -0 .95-0 .81-0 .27 




B~ 


L> 7T° 


T 




_ _„_i_n ^"-un n"-i-n 07 
u - ' °-0. 27-0. 06-0. 07 




B~ 


^ D K" 


A 


< 5.0 


~ 1 n +0. 06+0. 19+0.02 
U.19_o .06-0.16-0.02 




B~ 


-> D~r\ 


T 




n 07+0. 19+0.03+0.04 
-0.14-0.03-0.03 




B~ 


— > D~T]' 


T 




„ r)/) +0. 13+0.02+0.02 
U.^4_ _ 10 _o .02-0 .02 




B~ 


r\— o 

-> - 


T 




<-, ^i-t-Q 44-1 04-1 ft. 

1 ' -^-6.8-1.01-1.7 




B~ 


-» L> s K u 


A 




/r\ />o+0 39 + 90+0 07 \ -i r\ 


_2 


B~ 


-> D7v 


T 




-i /-i o+4 8+0 57+1 
lu -°-3.7-1.4-0.96 




R — 


. n — n' 


1 




r 71 +3.41+0.73+0.56 
di ' 1 -2. 28-0. 46-0. 53 




B° 


-» DV 


E 




(9.54l^l?; 3 ™) x 10- 


-2 


m 




C 




|-n ,0+1.4+1.6+0.24 \ |a- 
(, z -^°-1.0-0.98-0.23 ) x iu 


-2 


m 




c 




1 i n +0. 50+0. 25+0. 11 
x - lu -0. 41-0. 26-0. 10 




m 


-►ay 


c 




9 ro+1. 2+0.59+0.25 
z -°°-0. 94-0. 54-0. 24 




m 


-> D-7T + 


E 




n i q+0. 07+0. 01+0. 01 
U.l»_ . 06-0. 03-0. 01 








T 




1 91+0.58+0.10+0.12 
!- zl -0. 43-0. 13-0. 11 




m 


- ^7^+ 


T 




9Q 4+14.3+2.7+2.9 
^»-1_10.6-3.5-2.7 





decay constants, and some parameters of PQCD approach, are listed in the appendix [SJ With the formulae we list 
in Sec. UH we can get the amplitude A for each channel. The decay width for the — > £)? -.V mode is given by 



r = ^—(l-r*) J2 \^ X \ 2 (X = N,T) 



WlTUlB 

en,ev 

- T--— (l-r 2 )(|^| 2 + 2(|^| 2 + |^| 2 )). (52) 
lbirrriB 

For the other three kinds of decay modes B — ► D^P and B — > D( S )V, the decay width is given by 

-(l-r 2 M 2 - (53) 



WirrriB 

With the decay width at hand, the branching ratio is given by B7Z = Ttb- We take tr- = (1.674 x 10~ 12 )s/?i, 
t B o = (1.542 x 10~ 12 )s/?i, t B o = (1.466 x 10~ 12 )s/?i, and G F = 1.16639 x 10- 5 GcV -2 [13j. Our numerical results 



11 



TABLE II: Branching ratios of Br s ) — » DV decays calculated in PQCD approach together with experimental data[13| (unit: 
lO" 6 ). 



Modes 


Class Experiments 


Our results 


13 


n- + 
-> D s p T 


T 


< 600 


qc q+17.8+2.0+3.5 
OO.y_ 13 .2-1.9-3.3 


50 




t n 




(1 co+0. 63+0. 42+0. 16 \ in— 2 
(.I.DO^q 43_ .3i_o .15 ) x ll J 


Li 


n0 „0 
— > L> p 


C 




(I 44+1-8+1.2+0.34 \ ln -2 
V J -^-1.3-0.78-O.32 J x lu 


Li 


— > V K 


c 


< 11.0 


1 QO+0.95+0.75+0.18 
x - MZ -0. 77-0. 52-0. 17 


Li 




c 




(e.lSt^ili.n.s?) x 10- 2 


50 


— > U p 


T 




1 9^+0. 65+0. 07+0. 12 
1 - z ' -0.48-0.06-0.12 


D — 
15 


—> D p 


c 




(9.27^™) x lO" 2 


/ > — 


r>0 f*- 
D K 


i 




o n^+i-i+O-^+O- 20 

z - UJ -0. 76-0. 49-0. 19 


_/_.> 


u p 


1 




fl 7CI+0. 37+0. 05+0. 07 
u -' d -0. 28-0. 05-0. 07 


£5 


— > U i\ 


A 

A 




Oil + '~ , .Uo+U.Ul + U.Ul 

-0.04-0.02-0.01 


D — 
13 


^ 7~) — > < 


1 




r\ /?jy-|-U.oZ+U.Uo+U.UD 

U - D ' -0.25-0.04-0.06 


11 ~ 




1 




i n k i i aii n 

in y^+y.o+i.u+i.y 

iy -^-7.2-l.l-1.8 


11 ~ 
Li 


— > D s K 


A 

A 




> IA1(3IA1/1IAAK^ 1 

fn cn+ u -i b + u - i4 + u - Uj > v in -2 

l i U.OU_ . 15 _ .14-0 .05 J x iu 


11 ~ 

13 


— > i>> s W 


X 




■xr- o+8. 3+0. 94+1 .6 
rO.O_g .2-0.93-1.6 


11 ~ 
L) 


_^ D — 

U s 9 


A 

A 


1 Qfl 


1 Q + 0-UO+0.U3+0.U2 

u - lo -0. 04-0. 02-0. 02 


m 


-» 5V 


E 




fl 1 Q+0.05+0.04+0.02 
u - iy -0. 05-0. 03-0. 02 






C 




CO fin +0 - 32+0 - 27+0 06 ~1 X 10~ 2 

^U.DU_q 22— 0.17-0.06 ) X J-U 


m 




E 




n 1 fi+0. 05+0. 02+0. 01 
u - 1D -0. 05-0. 03-0. 01 


B° 




C 




1 OQ+1.0+0.71+0.18 
1 - o3 -0. 69-0. 45-0. 18 






E 




n 07+O.II+0. 06+0. 04 
u '°' -0.10-0.07-0.03 


m 


-> 


T 




1 ^9+0.63+0.08+0.14 
x -^ z -0.48-0.12-0.13 


m 




T 




00 1 +15.4+1.63+3.06 
•J - 1 -!!. 6-2. 65-3. 07 



are listed in Table III ITU IIIII and IIVI The first error is from the hadronic parameters of B^ meson wave functions(thc 
decay constants and the shape parameters): f B = (0.19 ± 0.025)GcV, f B „ = (0.24 ± 0.03)GeV, oj b = (0.40 ± 0.05)GeV 
for B meson, and cot = (0.50 ± 0.05)GeV for B s meson. The second error arises from the choice of the hard scales, 
which vary from 0.75t to 1.25i, and the uncertainty of Aqcd = (0.25 ± 0.05)GeV. This uncertainty characterize the 
size of the next-to-leading order QCD corrections, which is shown reasonable in these tables. The third error comes 
from the uncertainties of the CKM matrix elements: 

\V ub \ = 0.003 59 ± 0.000 16, \V cd \ = 0.225 6 ± 0.001 0, 

\V CS \ = 0.973 34 ±0.000 23, 7 = (77i|°)°- (54) 

It is easy to see that the most important theoretical uncertainty comes from the non-perturbative hadronic parameters, 
which can be improved later by the experiments. 

At quark level, the decay channels related in this paper are all b — > u transitions. This type of decays are suppressed 
by the CKM matrix elements, especially for the decays without a strange quark in the final states. This is one reason 
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TABLE III: Branching ratios of B( s \ — > D*P decays calculated in PQCD approach together with experimental data [l3j](unit: 
1(T 6 ). 



Modes 


Class Experiments 


Our results 




60 
Li 


n*- + 


T 


60. ± 7.U 


/II 7+21-9+2.5+4.1 
' -15.8-2.3-3.8 




60 
Li 




E 




(ft QC+0.10+0.16+0.04\ ln - 
V u - JU -0. 09-0. 10-0. 03^ A lu 


2 


50 
Li 


r.*0 


C 




^• ±9 -l. 61-1. 16-0.39 ) x iu 


2 


60 


— > L> A 


c 




9 qr + 1.2+0.94+0.23 
Z.JO_ 93_ .61-0 .22 




60 


— > D rj 


c 




/ /1Q+3. 86+2. 04+0. 82 \ ia- 
\°-^°-2. 94-1. 86-0. 78,1 x lu 


2 


60 


7=1*0 / 

— > L) rj 


c 




(5.64±; : g?±i;|5±g : 5|) x 10- 


2 


50 

B 


7~l*— + 
— > Z? 7T^ 


T 




1 91 +0.67+0.08+0.12 
1 - zl -0. 49-0. 07-0. 12 




B 


ri*0 - 
— > D 7T 


(J 




(R 1Q+4-29+2.33+0.60-, 1 n - 
(,0.10_ 3 . 17 _! .84-0 .57./ x lu 


2 




. fS*0 7y — 






ft 71+0.54+0.54+0.07 
u - ' 1 -0. 36-0. 39-0. 07 




B 


7~1*— 


T 


< 3.0 


ri ep+0. 36+0. 05+0. 06 
U - DO -0. 25-0. 05-0. 06 






. 7-)*- fT° 


A 

A 


< y.u 


-1 ^^+0.44+0.04+0.14 
i - 4D -0. 31-0. 10-0. 14 






— > U ?7 


1 




1 A 1 A m 1 A fir; 

u -°°-0. 19-0.05-0.05 




Li 


— > LJ rj 


T 
1 




ri qr + 0. 17+0. 02+0. 03 
u " 3O -0. 13-0. 03-0. 03 




Li 


— > 7T 


1 




f)t) 17+ i^.Z+ 1 .4+^.cJ 

zz - ' -8.7-1.2-2.2 




L> 


, 7~)*— 7/-0 






, 1 i nn 1 a 0*7 1 a •. 

/fi Qtr + ^-00+0.97+0.bH \ v in- 

l°- y0 -l. 65-1. 41-0. 65 J X 1U 




Li 


^ B a rj 


1 




1 /I S3A 1 A C1 1 1 A 7 
Q Q/1 +4.SU + U.OO + U.O 

°- y ^-3. 56-0. 93-0. 83 




D — 
/ J 


— ► U s 7/ 


1 




-LO-U_4. 42-0. 79-1. 21 








E 




1 A AC 1 A AO 1 A A1 
P) 1 O + O.OO+O.OZ+U.01 

u - 1 ^-0. 04-0. 02-0. 01 




B° 




C 




(3.25^:^™) x IO" 2 




B° 




C 




1 ic+o. 57+0. 21+0. 11 

± - ±u -0. 47-0. 26-0. 11 




B° 


-» D*y 


c 




n /M + l. 53+0. 71+0. 33 
25-0.71-0.32 




B.° 


-»■ D*~TV + 


E 




n 97+0.09+0.04+0.03 
u ' z ' -0.09-0.05-0.02 




B a ° 


-v D*~ K + 


T 




1 QO+0.66+0.11+0.13 
i -°°-0. 50-0. 16-0. 13 




B s ° 


-» d;- a-+ 


T 




oc 9+18.7+3.4+3.8 
JO - z -13.1-4.4-3.6 





why most of the decays have small branching ratios with order 10~ 6 or 1CP 7 . Another reason is the absence of the 
color allowed emission diagrams with a light meson emitted. 

As stated in previous section, all these decays do not have contributions from the penguin operators. For the 
tree operator induced decays, we have only four type of topology diagrams contributed: the color allowed diagrams 
(T), the color suppressed diagrams (C), the W annihilation decays (A) and the W exchange decays (E). All the 
decays are thus classified in the tables. From the numerical results, we can see that the pure annihilation type ("W" 
or "E" ) decay branching ratios are suppressed comparing with the "T" or "C" emission diagrams dominant decay 
channels. Remember that the "T" and "C" emission diagrams dominant decay channels may also have "W" and "E" 
annihilation type contributions, although they are suppressed relatively Within each category of decays, the large 
differences between channels are due to the Cabbibo suppression factor V c d/V cs . 

Usually the nonfactorizablc emission diagrams are suppressed comparing the factorizable diagrams. The two non- 
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TABLE IV: Predicted branching ratios of Br 3 \ — » D*V decays(unit: 10 6 ) and the percentage of transverse polarizations TZt 
together with the experimental data[l3|. 



Modes 


Class Experimental BRs 


PQCD BRs 




50 

13 


n*— + 

-» D a p + 


T 




(•q n+33.9+4.49+6.65 
DO. 4_ 25 2 _5. 24-6. 34 


(33±g : t±i:*)% 


B 


-> D*~K* + 


E 




CI Q1 +0.44+0.43+0.19 \ 1f) -2 
I !' al -0.63— 0.50-0.18 I x iu 


(40l 6 ; 8 i 6 s 9 °)% 


50 
B 


r>*0 
— > D p 


C 




fl Q1 +0.13+0.04+0.03 
U.Ol_ . 11-0. 05-0. 03 




13 


— > L) K 


c 


< 40 


1 ^ c+6. 20+1. 76+1. 32 
10 - J -4. 80-2. 16-1.25 


(81±g;l±|;|)% 


B 


ri*0 
— » £/ O) 


c 




O O/j +0- 11+0.27+0.02 
u - z *-0. 09-0. 04-0. 02 


(72l?;^: 4 5 )% 


50 
B 


n*— + 


T 




9 OQ+113+0.18+0.22 
z - z3 -0. 86-0. 17-0. 21 


(34tS:?tl:S)% 


B 


n*0 „— 


C 




1 in+ ' 42 +°- 15 +°l 1 
i - lu -0. 34-0. 42-0. 10 


(88l 2 ;;^; 9 3 )% 


B~ 


fS*0 t^* — 

— » L> K 


c 




OG. 1+9-93+3.68+2.54 
^O- 1 - 7. 90-4.88-2. 43 


(88±i;t±l;S)% 


B 




1 

1 




i en+ ' 70 + ' 11 + ai5 

±.OU_q 62 _ i5_o. 14 


(42 +l.l + 0.5 5)% 


B 


v n*- S-'O 


A 




i n i^o i n i n 
r\ rjr+U.bo + U.^O + U.^^ 

z - zo -0. 55-0. 31 -0.21 


(87±g;g±S:S)% 


B 


— > JJ LO 


1 




i ni +0.52+0.07+0.10 

J -- UJ --o.38-o.o6-o.o9 


(25±i;|±tJ)% 


u — 
B 


. n*- «o 


1 




Q + 18. 2 + 2.41 + J. oy 
13.6-2.69-3.42 


(33±g:i±i:«)% 


TJ — 

B 


. n*- z^*o 


A 

A 




ri -i -i +0.04+0.02+0.01 
u - 11 -0.02-0. 03-0. 01 


(89l 2 ;^°-i 8 )% 


B 


— > D s w 


1 




q-i q+15.7+2.09+3.11 
oii J -11.8-2.42-2.97 


(33l°;^i; 4 6 )% 


B~~ 




A 


< 12.0 


q fiQ+0. 86+0. 55+0. 30 
°- u9 -0. 73-0. 77-0. 29 


(89t°;<»t°;°)% 


B° s 




E 




n i 7 +0.05+0.04+0.02 
U.l 1 -0.04-0.03-0.02 


(38±|:2±^ 6 )% 


m 




C 




n aq+0. 21+0. 09+0. 05 
"•^-0.16-0.08-0.05 


(84 +0.99_+3 s 9 )% 


B° s 




E 




n 1 c+0. 04+0. 02+0. 02 
u - ld -0. 04-0. 02-0. 01 


(37ltii 6 3 8 3 )% 


m 




C 




i i 7 +4.89+1.88+1.15 
LL - ' -4.08-2.37-1.09 


(78±J;1±|;|)% 


B° s 


-> D*-p+ 


E 




n nx+0. 10+0. 06+0. 03 
U-O<±_ 08-0. 07-0. 03 


(39l 2 :Il^i 4 )% 


m 


-<■ D*~ K* + 


T 




9 QC + 106+0.18+0.23 
Z - JO -0. 81-0. 23-0. 22 


(32t° ™ 6 )% 


m 


-> Drr+ 


T 




P.A Q+29.6+5.57+6.32 
D1.O_22.4_6.90-6.03 


(32l^:i: 8 5 )% 



factorizablc diagrams in FigQJc,d) give nearly canceled contributions if the emitted meson is a light meson. However, 
it's not the situation here when the D meson is emitted, c quark and the light quark are very different in the emitted 
D meson. As a result, the nonfactorizable emission diagrams also give non- negligible contributions. For example 
those channels with the £)(*)° meson in the final state are color suppressed. The Wilson coefficients for factorizable 
contribution £j„ t and nonfactorizable contribution Mint are ai = C\ + Ci/3 and C2/3, respectively. Since ai ~ C%/3, 
the £mt(a2) and .Mm^C^/^) give similar contributions. For these color suppressed modes, one can find that the 
annihilation diagrams can also give relatively large contributions. Our numerical results indicate that sometimes the 
annihilation diagrams' contributions do have the same order of magnitude as the emission diagrams. We also find 
that the twist-3 distribution amplitudes play an important role, especially in the annihilation diagrams. 

When the charged meson in the final state, the emission diagrams are the color favored ones, with the Wilson 

coefficients a± = Ci + Ci/3 and C\/3, for the factorizablc diagrams and nonfactorizable diagrams, respectively. In 
this situation the nonfactorizable diagrams Mint(C\/2>) are highly suppressed by the Wilson coefficient, comparing 
with the factorizable diagrams £ini(&i), £mt( a i) >> -Mint (Ci/3). This means that the dominant amplitudes with 
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£mt( a i) are nearly proportional to the product of DW meson decay constant and a B to light meson form factor. 
This type of B — > DP decays have a little smaller branching ratios than those corresponding B — * DV decays, since 
the form factors of B — > V are a little larger. However, in the color suppressed modes, £mt(ci2) ~ M.int(C2/2>), or 
pure annihilation type decays, the above conclusion is not satisfied. Because the large Mint{C-2/^) or the annihilation 
diagrams will bring more complicated situations. So the branching ratios of B — ► DV are not definitely bigger than 
those of B — > DP. From the Table HI and HT1 we can see that some branching ratios of B — > DP are bigger than those 
of .B — > DV, which is different from b — > cuq decays. 

For the B^ — > Dj^V decays, we also calculate the ratio IZt of the transverse polarization in the branching ratios, 
which is given by 

2(\A°\* + \A>>\ 2 ) 

T |^ Ar | 2 + 2(|^ s | 2 + |^P| 2 )- [ ' 

From Eq. (|44p -Eq. (|5ip . we can find that the transversely polarized contributions of the emission diagrams are sup- 
pressed by the factor r, and those of the W exchange diagrams are suppressed by the factor rry . That is the reason 
in category "T" and "E" decays, we have relatively small transverse polarization factions. 

For color suppressed emission diagrams (C), the factorizable contribution ^((a^) and non-factorizablc contribution 
M-mtiCil'S) are at the same order magnitude, none of which can give dominant contributions. The two diagrams 
of Mint (longitudinal) cancel with each other and the two diagrams of Mj nt and M. p int (transverse) strengthen with 
each other. So we can expect large transverse polarized contribution of the branching ratio. From table HV| we can 
see that all the "C" type decays have transverse polarization around 80%. 

For the W exchange type decays, the factorizable annihilation contributions arc suppressed by the Wilson coeffi- 
cients, thus the dominant contribution is from nonfactorizablc annihilation diagrams. One can find that the ratios 
of the transverse polarizations for the W exchange diagrams are around 40%. For the W annihilation type decays 
(A), the factorizable diagrams dominate the branching ratios due to the large Wilson coefficients. However, the two 
factorizable annihilation diagrams strongly cancel with each other in the longitudinally polarized case, while they 
strengthen with each other in transversely polarized cases. In addition, there is also cancelations between the fac- 
torizable and nonfactorizablc contributions for the longitudinal polarizations. Therefore the transverse polarizations 
take a far more larger ratio in the branching ratios, which can be as large as nearly 90%. 

The transverse polarization ratio do not depend on the variation of the CKM factors, since these kinds of overall 
factors cancel in the ratio. The uncertainty shown in table HVl are from the hadronic uncertainty and factorization 
scale. Although these uncertainties are small, it does not mean that the polarization ratio is stable. In fact, it is quite 
sensitive to the hadronic wave function shape of the final state meson [lj and the power corrections. 

B~ — ► D°(D°)K~ decays can be used to measure 7 angle, see Ref. [10|, where the ratio r = j^[g~Z^o^-]j is an 
important quantity. With the amplitudes we obtain, this ratio is given as 



\A(B — > DK )| _ nQ9+ o.oi2+o.oo3 
\A(B~ D°K~y ' 



n nQ9 +0.012+0.003 ( rf.\ 
u - UJZ -0. 003-0. 003' V JL V 



The first error comes from the choice of the hard scales, and the second error comes from the CKM matrix elements. 
The uncertainty due to the hadronic parameters are canceled in the ratio, thus the calculation of the ratio is more 
precise and stable than the individual channels. With the formulas in Ref. 10(, 7 can be measured experimentally 
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through 



2 r (b+ -> d 1 k+) + r (b- -> DiK-)} 9 r - ; — " — — , . 

Ri = — — - - = 1 + r 2 ± J 4r 2 cos 2 7 - A 2 cot 2 7, (57) 

r(£+ ^D°if+) + r(B- -> D°K-) V ' " v y 

where i = 1,2, = (I? ± D°)/y/2 corresponds to the two CP eigenstates, i?i is defined as two charge-averaged 
ratios for two CP eigenstates, and A = A2 — Ai with 

A = T[B+ - DtK+) - T(B ~ - DlK ) (58) 

1 T(B+ D°K+) + T(B- -> D°K~) ' 1 ' 



In Ref. 22j, authors give similar diagrams on the relation of Eq. ([ST)) . One can find that the sensitivity of 7 to the 
other quantities increases as the r decreases. When r ~ 0.1 the extracted 7 will be very sensitive to A. Thus our 
value of r may be too small for the current experiments. 

IV. SUMMARY 

In this paper, we investigate B{B S ) — > D^P, D^V, D* S ^P, D* S ^V decays under the framework of perturbative QCD 
approach. We analyze the contributions of different diagrams in the leading order approximation of the mjj /rag 
expansion. It is found that the nonfactorizablc emission and annihilation diagram are also possible to give a large 
contribution. However, the emission contributions are still dominant in the branching ratios. All the branching 
ratios referred are calculated and the ratios of the transversely polarized contributions in the B(B S ) — ► DJ^V are 
also estimated. We find that the transversely polarized contributions, which mainly come from the nonfactorizablc 
emission diagrams and annihilation type diagrams, are very large. In some channels, they are even dominant. The 
branching ratios of this kind of decays are around 10~ 6 and 10 -7 , which means the method of extraction CKM angle 
7 is not effective in experiments. 
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APPENDIX A: WAVE FUNCTIONS AND DECAY CONSTANTS 
1. Wave functions of Br s ) mesons 



The -B( s ) meson wave function are decomposed into the following Lorentz structures 

Z e ifcl - z (0|6 Q (0)d /3 (z)|B (s) (P 1 )) 



15] 



(2nY 

= -^|( / p 1 + M Bw ) 7 6[0fl w (ft 1 )- J ^^B w (*i)]} • (Al) 

Here 4>B (s) {k\) and </>B (s) (fci) are the corresponding leading twist distribution amplitudes, and <^s (s) (fci) contributes 
little, so we neglect it. The final expression becomes 

$s (s) = -y==(/Pi + MB {s) )lh4>B {B) (ki). (A2) 
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The first determination of B- meson wave function was done in [16] . In our work the distribution amplitude in the 
b-space is a little differently: 



4>B (s) {x,b) = N Bia) x 2 (l - x) 2 exp 



2 l w fc j 2 



(A3) 



Here b is the conjugate space coordinate of kij_. Nb ( , ) is the normalization constant, which is determined by the 
normalization condition: 



/ dxd>B, ,(x,b = 0) 
Jo {> 2^J2N C 



(A4) 



2. Wave functions and decay constants of light pseudoscalar mesons 



The decay constant of the pseudoscalar meson is defined as: 

(0|?i7„75&|iW) =i/piV (A5) 
The Lorentz structure of light cone distribution amplitudes (for out-going state) for light pseudoscalar mesons is: 

(P(P 3 )ki Q (0)g 2/3 (*)|0> (A6) 



dxe 



ixP^-z 



[ 75 /P<f) A {x) + l5 m <p p {x) +m 7 5 (/> jh - 1)4> T (x)] 



where u is the light cone direction along which the light pseudoscalar meson's momentum is defined, and n is just 
opposite to it. The chiral scale parameter rriQ is defined as mo 



Mr 



m qi +m q2 ' 

The distribution amplitudes are expanded by the Gcgcnbauer polynomials and their expressions are 



4>T( X ) 



V2Nc 
fp 
2^/2N c 

fp 



3 x(l -x)[l + afcl /2 (t) + 4cl /2 {t) + a£cl /2 (t) 



c\' 2 {t)+4cy 2 {t) 



2y>2N c 

with t = 2x — 1 . The coefficients of the Gcgcnbauer polynomials are 



(A7) 
(A8) 
(A9) 



17] 



<2tt 



<2tt 



'37T 



0.44, 0^ = 0.25, 0^ = 0.17, 



l 2K 



0.2, 



The decay constants are 



0.43 , = 0.09 , a\ K = 0.24 , a\ K = -0.11 
0.55 , al K = 0.35 . 



f„ = 131MeV , Ik = 160MeV. 



(A10) 



(All) 



3. Wave functions and decay constants of light vector mesons 

The decay constants for the vector mesons are defined by 



(0|gi7„d^(P3,e)) - fvm v e^ {0\qia^q 2 \V(P 3 , e)) = ifi{e^P 3v - e„P 3/1 ). 



(A12) 
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Up to twist-3 the distribution amplitudes are 



(V(P 3 ,e* L )\q la (0)q 2 p(z)\0) = - 



V2N^ Jo 



dxe 



(V(P s ,e* T )\q la (0)q 2 p(z)\0) = --== J dxe***" [M v &<jfy{x)+ & /P 3 ^W 

+M v ie^ pal5 ^e* T 1 ' n p v^ v (x)} a/3 , 



(A13) 



where x is the momentum fraction of the q 2 quark. Contrary to the pseudoscalar case, here n defines the light cone 
direction along which the momentum of light meson is taken and v is just the opposite light cone direction. The 
twist- 2 distribution amplitudes of vector mesons are defined as 



M*) = rjf^ 1 - *) [ x + a " °l' 2 (*) + 4c 2 3/2 (t) 



3/v 



x(l-x) 1 + a^Cl /2 (t) + a^(J 2 r \t) 



-L/-.3/2, 



and the corresponding values of the Gcgcnbaucr moments are l^f 



*2p 



aL = 0.15 ± 0.07 , a\ K , = 0.03 ± 0.02 , a) 2 l K , = 0.11 ± 0.09 



a^ p = a 2ui = 0.14 ±0.06 , a^ K , = 0.04 ±0.03 , a 2K » = 0.10 ±0.08 , a 2<t> = 0.14 ±0.07 



l 2(j> 
2 2<j> 



0.18 ±0.08 



For the other distribution amplitudes, we use the asymptotic form 

3/£ 4.2 



(,<■) 



3/' 



^) = f^(l + i 2 )= ^(*) = ^|(-t) 



The decay constants are listed below: 



f p = 209 ± 2MeV, f K . = 217 ± 5MeV, f u = 195 ± 3MeV, / = 231 ± 4MeV, 
/J = 165 ± 9McV, /£. = 185 ± lOMcV, /J = 151 ± 9MeV, /J = 186 ± 9MeV. 



(A14) 



(A15) 



(A16) 



(A17) 



4. Wave function of Z)'*' meson 



Up to twist-3 accuracy the two-particle light-cone distribution amplitudes of meson are defined as 



(D(P 2 )\q a (z)cp(0)\0) = 



(D*(P 2 ) \q a (z)cf,(0)\0) 



dxe lx ^ z b 5 {/p2 + M)4> D {x,b)] afj 



i= f dxe^ z [/k L (fr + M D .)(&.{x,b)+ M/P2 + M D *)<j>l*(x,b)] (A18) 

JVC JO 



with 



dx(f>D (x, 0) 



D 



dxc/) D , (x, 0) 



D- 



,o 2V2N C 'J 2^2N C ' J 

as the normalization conditions. In the heavy quark limit we have 



dx(f>J), (x, 0) 



Jl_ 

2V2W: 



, T r m c + m d T m c + m d 

Jd* - JD* 



M D * 



Id* - fo* ■ 



0(A/M D .). 



(A19) 



(A20) 
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Thus we will use /q, = /d» in our calculation. The models for the distribution amplitude for D meson we used in 
this paper is 

C Gen) (x,b) = -±=f D 6x(l-x)[l+C D (l-2x)]exp[^^}. (A21) 
The masses of meson we use are 

m D = 1.869GeV, m D - = 1.968GeV, 

m D . = 2.010GeV, m D ,- = 2.112GeV. (A22) 

We take f D = 207McV [3, C D = 0.5, u = O.lGcV for D meson, f Ds = 241MevQ, C D = 0.4, u = 0.3GcV for D s 
meson. We use the relations between fr> and fn* derived from HQET [2^: 



Id* = \l—f D , f Dr = J-^f D7 - (A23) 

APPENDIX B: PQCD FUNCTIONS 

The pQCD functions appear in (7))-([TU]) and |[14" )) -([5"1"] ) are listed as 

h i (x 1 ,x 2l bi,b 2 ) = K (^xix 2 m B bi) 

x [6(bi - b 2 )K {y/x2~m B bi) I (y^msh) 

+6{b 2 -b 1 )K a {^m B b2)h{^m B b 1 )] , (Bl) 
h a {x 2 ,x 3 ,b 2 ,b 3 ) = [i-J (y/x 2 x 3 m B b 2 ) 

b 2 - b 3 )H^ (y/x3~m B b 2 ) J (y/x^m B b 3 ) 
+6(b 3 - b 2 )H { 1} (^msh) Jo (Vxlm B b 2 )\ , (B2) 

where H^(z) = J (z) + iY (z). And 
,0') 



K> = [0(bi - b 2 )K Q (Dm B bi)I (Dm B b 2 ) + 9{b 2 - h)Ko (Dm B b 2 ) I (Dmsh)] 
Ko(Djin B b 2 ) for D] > 



H^\J\Df\m B b 2 ) forD?<0 



2 



h U) - i- 



with the variables 



(B3) 



0(b! - b 2 )H^ (FrriBh) J (Fm B b 2 ) + 9{b 2 - b^H^ (Fm B b 2 ) J Q (Fmsh) 

KoiFjmBh) for F 2 > 
ifH^{yJ^f\m B bi) forFf <0 ' 

D 2 = x 1 (l-x 3 )(l-r 2 ), 

D\ = (a; 1 -x 2 )(l-X3)(l-r 2 ) , 

D\ = (xi + x 2 )r 2 - (1 - Xl - x 2 )(l - x 3 )(l - r 2 ) , 

F 2 = x 2 (l - x 3 )(l - r 2 ) , 

F 2 = 1 - (I - x 2 )(l - Xl - (1 - x 3 )(l - r 2 )) , 

F 2 = x 2 ( Xl - (1 - x 3 )(l - r 2 )) . (B5) 
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The hard scales are determined by 



tW = 



max 



( v /(l-x 3 )(l-r 2 )m B , I/61, 1/63) , tf ] = max(v/xi(l -r 2 )m s , l/b u I/63) 



max 



( v/a; 2 (l - r 2 )m B , l/6 2 , 1/63) , = max( ^(1 - x 3 )(l - r 2 )m B , l/b 2 , 1/63) 



t)>> = max(Dm B , ^/|£>||m B , 1/62) , = max(Fm B , ^\Fj\m B , l/h, 1/b 
Jet function appears in the factorization formulae is 

*'M° 2 Tm! + ;W -*r. 

Vtt r(i + c) 

with c = 0.5 in this paper. 

The expressions for the functions Ei(t), with I = i,a, d, /, are 



(B6) 



(B7) 



= a s (t)exp[-S B (t)-S 3 (t)}, 

E a {i) = a s (t)exp[-S D (t)-S 3 (t)}, 

E d (t) = a s {t)cxp[~S(t)\ b3=hl ] , 

E f (t) = a s (t)exp[-S(t)\ b3=b2 ] , 



(B8) 



where the Sudakov exponent S = Sb + So + S3, with S3 as the Sudakov exponent of the light meson. And the 
Sj(xi)(j — B,D or 3) functions in Sudakov form factors are 

H 



m B 



S B (t) = s [ si^p&i ) +2 



dfi 

1/61 Z 2 



7g(a s (»), 



V2 



-^7«(a s (M)) ) 
1/62 A* 



S 3 (t) = s [x 3 -^,b 3 )+s (l-aa) -£,b. 



V2 



m B 



V2' 



^7g(a s (^)), 
1/63 M 



(B9) 
(BIO) 
(BH) 



with the quark anomalous dimension 7 9 = —a s /ir. The explicit form for the function s(Q,b) is 

Ad 



s(Q,b) 



gin 



2/?i 
A (1 )/? 2 



AC 1 ) / f N A( 2 ) 



2/3i 



6 

ln(2g) + 1 ln(26) + 1 



m \b 



1 



8/3? 



A( 2 ) Ad) /e 2 ^" 1 
~ In 

4/3 2 4/3i 



ln 2 (2<7) -ln 2 (2&) 



In 



(B12) 



where the variables are defined by 



q = ln[Q/(V2A)], b = ln[l/(6A)], 



(B13) 



and the coefficients A^ and are 



33 - 2n f 153 - 19n f 

Pi — — 1 P2 = - 



12 

-2 
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Ad) = l, A(3) = ^-^-H n/ + Vln(V), 
3 9 3 27 J 3 M y 2 h 



(B14) 



n/ is the number of the quark flavors and 7_e is the Euler constant. We will use the one-loop running coupling 
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constant, i.e. we pick up the four terms in the first line of the expression for the function s(Q, b). 
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